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Abstract 

Using data collected by the CLEO III detector at CESR, we report on the first observation of the 
decays — > A7r"'"7r"'"7r~7r'' and A+ Awtt"*". The branching fractions are measured relative to the 
pK'-K-^ mode; we find B(A+ ^ A7r+7r+7r-7rO)/fi(A+ ^ pK-TT+) = 0.36±0.09(stat.)±0.09(syst.) 
and B{A+ AuTr+)/B{A+ pK-Tr+) = 0.24±0.06(stat.)±0.06(syst.). We also observe a small 
excess of events in A^ Ar/vr"*", r/ — ^ 7r"'"7r~7r'^ from which we find -B(A^ A7]tt~^)/B{A'^ — > 
pK~TT~^) = 0.41ib0.17=b0.10 with a corresponding limit of <0.65% at 90% confidence level. This 
rate is consistent with previous measurements of A^ ^ A?77r+, r] — > 77. We find that the combined 
branching fractions for A+ Alott'^ and A+ Arjir^ comprise (87±13)% of the total A+ 
A7r"'"7r'^7r~7r° rate. We set an upper limit on B{Af — A7r+7r'^7r"'7r°)/iB(A+ — pK''Tr'^) for non-o; 
and non-?7 modes of 0.13 at 90% confidence level. 

PACS numbers: 13.30.Eg, 14.20.Lq 
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Over the last decade large data samples have been used to measure charm meson and 
charm baryon branching fractions to exclusive final states. Exclusive final states for 
A(n)7r decays, with n=l,2,3 have been measured but modes with four (or more) pions 
have yet to be observed. The large rates for K'^Tr^n^n^ir^ (5.4tf;4)%) and 

D° /^^Tr+TT+Tr^Tr" ((4.0±0.4)%) suggest non-negligible rates for A+ A7r+7r+7r"7r°. In 
this paper, we report the first measurement of the branching fractions for A^ —>■ A7r+7r+7r^7r° 
and A+ —>■ Acutt"'", uj —>■ tt'^tt'ti^. There is also some indication that a portion of the 
A7r"'"7r+7r^7r° occurs through A^ Ar^vr"*", rj — > 7r+7r^7r°. We find that the resonant u and rj 
sub-modes nearly account for most (if not all) of the A^ Ati^ti^ti^h^ rate. Throughout 
this paper, charged conjugate states are implied unless otherwise specified. 

The data for this analysis were collected using the CLEO III detector located at the 
Cornell Electron Storage Ring (CESR). The data sample consists of 2.7 fb~^ on or near the 
T(4S), 1.4 fb-i on or near the T(3S), and 0.02 fb'^ on or near the T(1S). The CLEO III 
detector is an upgrade of the CLEO II/II.V detector |^ which includes a four-layer double- 
sided silicon vertex detector, a tracking drift chamber [|] consisting of 16 inner axial layers 
and 31 outer small-angle stereo layers, and a ring-imaging Cerenkov detector (RICH) The 
RICH is located outside the tracking chambers and uses LiF radiators to produce Cherenkov 
photons from incident charged tracks. The Cerenkov cone expands in a 20-cm thick N2- 
filled expansion volume and are then photo-converted in a gas mixture of Methane and 
TEA and detected by a multi-wire chamber with cathode pad readout. The electromagnetic 
calorimeter consists of 7800 Csl crystals which are used to measure the position and energy 
of photons and to aid in lepton identification. These detector components are immersed in a 
1.5 Tesla solenoidal magnetic field that is uniform to within 0.2% over the tracking volume. 
Outside the solenoid there is a muon system covering 85% of the solid angle. 

Signal candidates are formed using charged tracks and photons. We only use charged 
tracks which have a high quality track fit and have | cos^| < 0.925, where 6 is the polar 
angle of the track with respect to the e~ direction. Charged particles (except those used in 
forming A candidates) are required to have a distance of closest approach (to the interaction 
point) in the r — (f) plane of 5 mm and 5 cm along the z axis (the direction of the e~ 
beam, apart from the small crossing angle of the e^e~ beams). Photons used in forming 7r° 
candidates must be isolated from charged tracks and have a shower shape that is consistent 
with the shape of an electromagnetic shower. The isolation criterion requires that no charged 
track points to any of the crystals which are used to make up the electromagnetic shower. 

Particle identification is performed using information on the specific ionization (dE/dx) 
measured in the drift chamber and particle likelihoods obtained from the RICH. For dE/dx, 
we quantify the consistency of the measured dE/dx with p, K and vr hypotheses by defining 
a significance. Si = {{dE / dx)f'^ — {dE / dx)i) / ai{dE / dx), where {dE/dx)f'^ and {dE/dx)i are 
the measured ionization (corrected for path length) and expected (momentum-dependent) 
mean value of dE/dx for particle type i (i = p,K,7i) respectively, and ai{dE/dx) is the 
resolution on the measured dE/dx. For the RICH, we use the measured track momentum 
and a given particle hypothesis to search for Cerenkov photons whose measured positions 
are consistent (within 3 standard deviations) with expected positions. For each particle 
hypothesis, z, we compute a likelihood of the form: 

u= xHGiefieT") + B]. 

Here, G is a Gaussian-like function which expresses the probability of observing the j^^ 
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photon at an angle 9° with respect to the expected Cerenkov angle 9^^^ for particle type i. 
The second term, B, is the background probability distribution which is a flat function in 
Qexp _Qobs^ Using the likelihoods, we compute differences in negative- log likelihoods between 
two particle-type hypotheses: proton-kaon, proton-pion, and kaon-pion (i.e., C{K — vr) = 
— 2 log(L;^/L7r)). Defined in this manner, true kaons will peak at negative values of C and 
true pions will peak at positive values of C 

Proton candidates are identified by requiring either \Sp\ < 3, 15*^1 > 3 and \Sk\ > 3, 
or at least three reconstructed Cerenkov photons with C{p — tt) < —5 and C{p — K) < 
0. Kaon candidates are only required to have either \Sk\ < 3 and \S-„\ > 3, or at least 
three reconstructed Cerenkov photons with C{K — vr) < 0. Because pions are the majority 
particle in hadronic events, only loose particle identification requirements are imposed. Pion 
candidates are required to have 15,^1 < 4 or at least three reconstructed Cerenkov photons 
with C{K ~ n) > —5 and C{p — tt) > —5. By construction, the list of pion, kaon and proton 
candidates are not exclusive of each other, but in forming a given candidate, all daughter 
tracks must be distinct from one another. 

Candidate A baryons are constructed by pairing oppositely charged proton and pion 
candidates and performing a constrained vertex fit. To suppress combinatorial background 
we require A candidates to have a decay length of at least 1 cm. The invariant mass 
distribution of A candidates is shown in Fig. 0(a). The distribution is fit using two Gaussians 
(with the same mean) to describe the signal and a first order polynomial to describe the 
background. The relative fraction of events in the Gaussian components is consistent with 
simulation. The signal (sideband) region includes candidates which have a mass in the range 
from 1111-1121 MeV/c2 (1100-1105 or 1130-1135 MeV/c^). The A sidebands are used to 
check for a possible selection bias. Candidate mesons are formed by combining pairs of 
photons and requiring that the pull is less than 2. The pull is defined as the difference in 
the 77 invariant mass from the known vr" mass, divided by the uncertainty in the 77 mass, 
where the uncertainty on the invariant mass is computed on a combination-by-combination 
basis. The pull distribution for the candidate tt^'s in the A+ ATr+Tr+Tr^Tr^ sample are 
shown in Fig. |l](b). Superimposed is the result of a fit to the sum of a Gaussian and a 
linear background. The fit has a mean of zero and a width of 1.1 units which is consistent 
with simulation. To reduce combinatorial background, we require that the energy of each 
photon is larger than 50 MeV and that the normalized energy asymmetry of the two photons 

1^1-^21/(^1 + ^2) < 0.8. 

Candidates for the decay A^ — A7r+7r+7r^7r° are formed by combining A candidates with 
one neutral and three charged pions. The combinatoric background is suppressed by re- 
quiring A^ candidates to have a total momentum in excess of 3.5 GeV/c. Kinematically, 
this removes the contribution of A+ baryons from B meson decay, leaving only the contri- 
bution from charm fragmentation. To improve the signal-to-noise ratio, we require A, 7r° 
and candidates to have momenta in excess of 500 MeV/c, 250 MeV/c, and 100 MeV/c 
respectively. 

Candidates for the A+ pK~7T~^ decay (reference mode) are formed by combining proton 
candidates with charged kaon and charged pion candidates. As with the A+ — > A7r"'"7r"'"7r~7r° 
candidates, we require that the A+ —>■ pK~7r^ candidates have total momentum larger than 
3.5 GeV/c. The charge of the kaon is required to be opposite to that of the proton. As the 
decay length of the A+ is too small to be observed in the CLEO detector, we require that 
all three tracks are consistent with coming from the e~^e~ interaction point. We also require 
that the proton have a total momentum larger than 500 MeV/c, while the kaon and pion 
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are required to have momenta larger than 300 MeV/c. 

Figure ^(a) shows the invariant mass distribution for A+ — A7r+7r+7r^7r° candidates for 
data after all selection criteria. A clear excess of events is observed at the A+ mass. The 
distribution is fit to the sum of a Gaussian whose width is fixed to the expected width from 
simulation (8.2 MeV/c^) and a linear background from which we extract a signal yield of 
49.5±11.9 events at a mass of 2287±1.3 (stat.) MeV/c^. Also shown in the figure (shaded 
histogram) are combinations where we pair the 7r+7r+7r~7r° with a A. No A+ signal should 
be present in this mode, and the data show no significant excess in this mass range. The 
data have also been separated into A+ and A+ in Fig. ^(b) and (c), respectively. They are 
fit using the same functional form as for the entire signal. The extracted yields are 21.3±8.1 
A^ and 28.3±8.5 A+ events. We therefore confirm that A^ and A+ appear in approximately 
equal numbers. 

We search for substructure in the A7r+7r+7r^7r° final state by forming tt^ti^ti^ combi- 
nations (2 per event) among the 4 pions using A7r+7r+7r^7r'^ combinations which have an 
invariant mass between 2270 and 2300 MeV/c^. The 7r+7r~7r'' invariant mass distributions 
are shown in Fig. ^ for the A+ signal region (2270-2300 MeV/c^) and the A+ sideband regions 
(2245-2260 and 2320-2335 MeV/c^). A strong peak is observed near = 782 MeV/c^ and 
a smaller excess at = 547 MeV/c^. The distribution is fit to the sum of two Gaussians, 
one for the u and one for the t], and a second-order polynomial for the background. The 
signal shapes used in the fit are determined from simulation. The fitted masses are 552±2 
MeV/c^ and 778±3 MeV/c^, consistent with the known masses of the rj and u mesons. The 
probability that the enhancement at the t] mass is caused by a background fluctuation is 
about 4%. The u signal is approximately 6 standard deviations above the expected fitted 
background. We also observe three narrow peaks between 825 and 900 MeV/c^. Each of 
these peaks correspond to ~1.5-2 standard deviations above the background prediction, de- 
pending on the assumed Gaussian width. Since there are no known particles at these masses, 
we assume they are statistical fluctuations. 

The signal yields in the uj [rj) modes are determined by requiring the 7r"'"7r~7r*^ invari- 
ant mass (either combination) to fall within ±30 MeV/c^ (±14 MeV/c^) of their known 
masses and fitting for a signal in the A7r"'"7r+7r~7r*^ invariant mass spectrum. To ensure non- 
overlapping sub-samples, we exclude A^ —>■ Atjtt^ candidates from the A+ — > Autt^ sample. 
The signal regions for the uj and 1] correspond to about ±3.5 times the Gaussian width 
obtained from simulations of these decay modes. The A7r+7r+7r~7r° invariant mass spectra 
for the u and rj modes are shown in Fig. ^(a) and (b), respectively. The shaded regions show 
the corresponding distributions when we select combinations which fall into the u sidebands 
(690-720 MeV/c2 or 843-873 MeV/c^) and r] sidebands (500-513 MeV/c^ or 580-593 MeV/c^) 
respectively. The distributions are fit to the sum of a Gaussian and a linear background. 
The widths of the Gaussians are constrained to 8.4 MeV/c^ for the u and 6.6 MeV/c^ for 
the 7]. These widths are determined from Monte Carlo studies of these processes, which 
includes the Breit-Wigner width of the u and detector resolution effects. The fitted yields 
are 31.9±7.1 A+ Aujtt^ candidates and 10.8±4.1 A^ Arjir^ candidates. Due to the 
limited statistics in the A^ Arjir^ mode, we also quote an upper limit of 16.8 events at 
the 90% confidence level. Thus these two modes account for 87% of the A+ —>■ A7r"'"7r+7r~7r° 
candidates. 

The contribution from the uj and rj are combined and shown in Fig. ^(a). The spectrum 
is fit to the sum of a Gaussian and a linear background shape. The width is constrained to 
the event-weighted average of the widths obtained from the Ajt — ^ Acutt"*" and A+ —>■ ArjTT^ 
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simulations (7.9 MeV/c^). The number of fitted candidates is 42.5±8.2, which agrees with 
the sum from the individual fits. Events which do not have a '7r"'"7r~7r'' with invariant mass 
in the t] or u region are shown in Fig. ^(b). We fit the spectrum assuming a Gaussian signal 
shape and a first order polynomial for the background shape. The mean of the Gaussian is 
fixed to the world-average value (2285 MeV/c^) and the width is constrained to the value 
found from a simulation of non-resonant — >■ A7r"'"7r"'"7r~7r'^ (9.8 MeV/c^). The fitted yield 
of S.Oly;^ events is consistent with the expected difference of 7.0 events as well as zero. This 
corresponds to an upper limit of 17.7 events at the 90% confidence level. 

The fitted Aujtt^ and Atjtt^ signals are corrected for possible contributions 

from non-resonant production by subtracting off the expected non-resonant contribution, as 
determined from simulation. The number of non-resonant events is taken to be 7.0 events 
as discussed above. This number must be corrected for the contribution from A+ —>■ Aujtt^, 
because some of the signal lies outside the u mass window. The rj width is negligible 
compared to the mass resolution, and therefore no such correction is necessary for A^ — »• 
Ari7r~^. From simulation, we find that (3±2)% of the u signal events fall outside the oo 
mass window. We therefore estimate that 6.0 of the 49.5 events are from sources other than 
A+ — >■ Au!7r~^ and A+ — >■ At^tt"*". We assume that these 6 events are from non- resonant sources, 
and using simulation we find that 16% (4%) are expected to have a vr'^Tr^Tr*^ combination 
in the u {rf) mass window. The corrected yields are therefore 31.0±7.1 and 10.6±4.1 for 
the A+ — i> AuoTi^ and A+ — > A?77r"'" modes. The 3% loss of uj events due to the mass 
window definition are accounted for in the A^ Acuvr^ efficiency. We take the systematic 
uncertainty on the non-resonant subtraction to be 100% of the correction. 

We have investigated the possibility that the uott system may come from higher mass 
resonance decays, such as A+ Ap'~^,p'~^ un^, which may contribute to this final state 
due to its large width. We compared the ujtt^ invariant mass distribution in data to a non- 
resonant tuTT"^ simulation and the expectation from p'^ —>■ ujtt^. No peaks are observed in 
data, and a Kolmogorov-Smirnov test of the shapes give a consistency of 80% with respect 
to a pure non-resonant Aujtt^ final state and a 10% consistency with pure Ap'"*". Thus, 
the data are more compatible with a pure non-resonant ujtt'^ contribution, although a p'~^ 
contribution cannot be ruled out. 

The relative branching fraction -B(A+ An^7i^7i^7T^)/B{A'^ —>■ pK^n^) is: 

B{A+ A7r+7r+7r-7r°) _ A^(A+ ^ A7r+7r+7r-7r°) e(A+ ^ pK-n+) 
B{A+ ^ pK~TT+) ~ N{A+ pK-TT+) e(A+ ^ ATr+Tr+Tr^vrO) ' 

where A^(A^ An^ 71^71^71^) and A^(A^ — > pK^TC^) are the number of reconstructed 
A+ — s> A7r+7r+7r^7r° and A^ pK^TC^ signal events and e(A+ — ATr+Tr+Tr^vr'^) and 
e(A^ pK~7c'^) are the relative efficiencies for reconstructing these modes respectively. 

The efficiencies for the Avr+vr'^vr^vr'^ and the pi^vr final states are obtained using a Monte 
Carlo simulation of these decays (see Table |). Events are simulated using the Lund/Jetset 
model for the production mechanism {e~^e~ qq) and are allowed to decay according to 
the known lifetime and decay branching fractions 0. The response of the detector to the 
particles in the generated events is simulated using a GEANT-based simulation 0| of the 
detector. The fully simulated events are then passed through the same analysis chain as data. 
For the pKi: final state, our simulation includes the resonant sub-modes pK* and A'^+i^ as 
well as a non-resonant pKir contribution 0. The efficiency for the A^ A7r+7r+7r^7r° final 
state is taken to be the branching fraction weighted efficiency for A^ Au;7r+, A+ Ar/vr"*", 
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and non-resonant A7r+7r"'"7r~7r°. Table | summarizes the inputs for the computation 

of the branching fractions. 

The branching fractions for A+ A7r+7r+7r^7r^, A+ Aujtt^ and A^ — > Arj-R^ relative 
to A+ pK^7i~^ are measured to be: 



^ ' ^ = 0.36 ±0.09 ±0.09 



BiA+ pK'7T+) 



B{A+ pK-n 



0.21 ±0.05 ±0.05 



B{A+ pK-n^ 



0.093 ± 0.038 ± 0.022 (< 0.15 at 90%CL), 



where the first uncertainty is statistical and the second is systematic. The systematic un- 
certainties are discussed below. 

These measurements are the first observations of A+ Ati^ti^ii^h^ modes, which are 
(or almost) entirely through A+ — > Aujtt'^ and A+ Arj-K^ . This is the first observation 
of the A+ — >■ AcuTT"*" decay mode. The A+ — Arjn^ decay has been previously observed 
in ?7 — >■ 77 decays |@], but this is the first indication in the rj — ^ tt^tt't^^ mode. Our 
measurement compares well with those results which yielded B{A^ Arj-K^) / B{A'^ — >■ 
pK--K+) = 0.35 ± 0.05(stat.)±0.06(syst.). 

As a systematic check of our procedure, we have also measured the relative branching 
fraction for A+ ^ A7r+ and find B{A+ A7i+)/B{A+ pK-7i+)=0.17±0.02{stat.). This 
is consistent with the current world average of 0.18±0.03 0. 

Significant systematic uncertainties arise from sources which affect the A7r+7r+7r^7r° mode 
but not the pK~'K^ mode (or vice versa). Systematic uncertainties include estimation of 
reconstruction efficiencies of tracks, A and 7r° reconstruction, particle identification selection, 
background estimation and signal estimation. 

Studies of data and simulation show that the charged particle tracking simulation re- 
produces data within 2%. A 2% systematic uncertainty per track introduces a 6% relative 
systematic uncertainty. The systematic uncertainty in the 7r° reconstruction efficiency is 
estimated by comparing the differential yield for T(45') — > vr^X in CLEO III data to CLEO 
II data. The difference is found to be 8% for our selection criteria. When combined with 
the 5% systematic uncertainty on the CLEO II 7r° efficiency, we assign a 10% uncertainty 
to the vr'' reconstruction efficiency. The systematic uncertainty in the A reconstruction effi- 
ciency is estimated by comparing our measured value for i3(A+ — ^ A7r"'")/i3(A+ pK~n~^) 
[0.17±0.02(stat.)] to the world average value (0.18±0.03). While these results are consistent 
with one another, we take the uncertainty to be the the quadrature sum of the difference 
in the two results (6%) and the relative uncertainty in the world average value (17%). The 
uncertainty in the A reconstruction efficiency is therefore estimated to be 18%. The sys- 
tematic uncertainty in the particle identification was investigated using a sample of protons 
from A — > pTT~ and kaons from D*~^ D^{^ K~7r~^)7r~^. The studies indicate that the 
systematic uncertainty is less than 3% for protons and less than 1% for kaons. The uncer- 
tainty in the background is estimated by using alternate background shapes derived from 
(a) the A7r"'"7r"''7r~7r° invariant mass distribution in data and (b) the A7r"'"7r"^7r~7r'' invariant 
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mass distribution using the A sidebands. We also varied the number of bins used in the fit 
to the background; these sources lead to an uncertainty of 6% on the signal yield. 

Additional uncertainty in the signal estimation may arise from an imperfect understand- 
ing of the signal resolution and efficiency. The resolution uncertainty is estimated by varying 
the widths by one standard deviation ± 1 MeV/c^) and refitting the data. The uncer- 
tainty in the signal yield is found to be 5%. For the efficiency uncertainty, we consider one 
source which only affects the full A+ A7r+7r"'"7r~7r° branching fraction measurement, and 
one which affects the resonant lj and rj contributions. For the full A^ A7r"'"7r"'"7r~7r'' sam- 
ple, we use a branching fraction weighted efficiency. To evaluate the systematic uncertainty, 
we shift the number of u events by ±7.6 events and the corresponding number of t] events 
by =f7.6 events and compute the corresponding changes in the A+ — > Att+tt+tt"??'' branching 
fraction. We find a shift of ±3% in the A^ — * A7r"^7r"'"7r~7r'' branching fraction. For the 
resonant modes, we also include a 3% systematic which corresponds to a 100% relative un- 
certainty on the subtraction for possible non-resonant contributions which have a tt+tt^tt^ 
invariant mass in either the u or rj mass region. A compilation of the sources of systematic 
uncertainty is presented in Table ||. The total systematic uncertainty is estimated to be 
13%. 

These measurements indicate that (87 ±13)% of the four-pion state is feed-down through 
resonant u and rj production. The branching fractions (after dividing through by the B{uj — > 
7r"'"'7r~7r°) and B{r] — >■ 7r"'"7r~7r°)) are shown in Table |ITT[ The common systematic uncertainty 
in the branching fraction for A+ — > pK~ 7r~^ has not been included in the last column in 
Table The A^ Arjir^ mode has been measured previously in the 77 — 77 mode [Q, 
and the result presented here is consistent with those findings. We find that the branching 
fraction for A+ Aw7r+ is (0.24±0.06±0.06) of i3(A+ pK--K+). For the r] mode, we find 
a rate relative to i3(A+ pK'7c+) of (0.41±0.17±0.10) with an upper limit of 0.65 at 90% 
CL. The non-uj , 1] signal constitutes S.Oty]/ events which corresponds to an upper limit of 
17.7 events at 90% confidence level. This corresponds to an upper limit on the branching 
fraction relative to B{A^ pK^n^) of 0.13 at 90% confidence level. 

The fact that the A7r+7r+7r~7r° rates are competitive with Avr^ ((0.9±0.3)%), Att+tt^ 
((3.6±1.3)%) and Att+tt'tt"'' ((3.3±1.0)%) is not surprising, as charm mesons exhibit similar 
branching fractions for D —>■ Kn{7r), where n = 1 — 4. In fact, it is interesting to note that 
the exclusive rate for D° K^un ((2.7±0.5)%) constitutes about 68% of the rate for 
— »• i^'~7r+7r+7r~7r° ((4.0±0.4)%). This is quite comparable to what we find here for the 
ratio of the exclusive Auir^, uj — > tt^h'tt^ rate to the total A7r"'"7r"'"7r~'7r° rate (59%). 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. M. Selen thanks the PFF program of the NSF and the 
Research Corporation, and A.H. Mahmood thanks the Texas Advanced Research Program. 
This work was supported by the National Science Foundation, and the U.S. Department of 
Energy. 
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TABLE I: Inputs for the computation of the branching fractions for A^ A7r"'"7r"'"7r~7r'', A^ — > 
Acutt"*" and A^ Arj-K^ . The label "tot" refers to the total A^ Att^tt^tt^tt^ sample for which 
the efficiency is computed as a weighted average efficiency of efficiencies as described in the text. 
The label "nr" refers to the simulated non-resonant A+ A-Tr+Tr+Tr^Tr^ sample. 

Mode Data Yield Eff. (%) g (7r+7r'7r")(%) 

pK^ 1959±57 20.8±0.9 - 

A7r+7r+7r-7r° (tot) 49.5±11.9 1.46±0.10 

AujTr+ 31.0±7.6 f.56±0.10 88.8±0.7 

A?77r+ 10.6±4.1 1.22±0.10 22.6±0.4 

ATT+TT+TT-TrO (ur) 7.9 1.48±0.10 



TABLE IL Systematic uncertainties in the measurement of the branching fractions for A^ 

A7r+7r+7r~7r'^, A+ Aa;7r+ and A+ Arjir'^. 

Source Value 

Background shape ±6% 

Signal Width ±6% 

Track Recon. Efficiency ±6% 

7r° Recon. Efficiency ±10% 

Proton ID efficiency ±3% 

Kaon ID efficiency ±1% 

A Recon. Efficiency ±18% 

Non-resonant subtraction (w, ij only) ±3% 
Signal Efficiency(A+ A7r+7r+7r-7r° only) ±3% 

Total ±24% 



TABLE III: Summary of measured A^ branching fractions. The uncertainties shown are statistical 
(first) and systematic (second). The qualifier "tot" refers to the total yield in the A7r+7r+7r~7r° 
final state and "nr" refers to non-a; and non-r? contributions. For the A+ Arjir'^ mode, we show 
both the central value and the 90% confidence upper limit (in parentheses). 





Mode 


;B(Mode)/;B(A+ ^pii:-7r+) 


;B(Mode) (%) 


A+- 


A7r+7r+7r~7r° (tot) 


0.36 ± 0.09 ± 0.09 


1.79 ±0.47 ±0.43 




A+ Au;7r+ 


0.24 ± 0.06 ± 0.06 


1.19 ±0.30 ±0.29 




A+ Ar/7r+ 


0.41 ±0.17 ±0.10 (< 0.65) 


2.06 ± 0.84 ± 0.49 (< 3.3 ) 


A+- 


-> A'!r'^TT'^'K~7r^ (nr) 


< 0.13 


< 0.65 
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FIG. 1: (a) Proton-pion invariant mass distribution after the proton identification requirement. 
The hghtly-shaded area corresponds to the A signal region and the darker shaded areas are the A 
sidebands, (b) The pull distribution for tt^ mesons which have momentum larger than 250 MeV/c 
and are in our A+ candidate sample. The superimposed curves represent fits to the distributions 
as described in the text. 
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FIG. 2: Invariant mass distributions for (a) both ATr+vr+Tr^vr'^ and Avr^Tr^Tr+vr'^ combinations, (b) 
Avr^vr^Tr^vr'^ combinations only, and (c) A7r~''7r"'"7r~7r'^ combinations only. In (a), the shaded distri- 
bution corresponds to candidates formed from the improper decay sequences (A^ — s- A7r''~7r'''7r~7r'^ 
and A^ — > Avr^vr^TT^Tr^). The superimposed curves are fits to a Gaussian signal function and a 
linear background as described in the text. 
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FIG. 3: Invariant mass distribution for the two tt+vt vr^ combinations in the Avr+Tr+Tr vr'^ 

decay for the signal region (histogram, 2270-2300 MeV/c^) and the A+ sidebands (shaded, 
2245-2260 and 2320-2335 MeV/c^). The superimposed curve represents a fit to the distribution as 
described in the text. 
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FIG. 4: Invariant mass distribution for A7r"'"7r+7r^7r'^ combinations from data where either of the 
two 7r"^-7r~7r'^ combinations have (a) a mass in the u) mass region (752 to 812 MeV/c^) or (b) in 
the r/ mass region (534 to 560 MeV/c^). The shaded regions show the corresponding distributions 
for the LO sidebands (690-720 MeV/c^ or 843-873 MeV/c^) and the rj sidebands (500-513 MeV/c^ 
or 580-593 MeV/c^). The superimposed curves are fits to a Gaussian signal function and a hnear 
background as described in the text. The fitted numbers of events are shown at the top of each 
figure. 
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FIG. 5: Invariant mass distribution for A7r"'"7r+7r^7r'^ combinations from data where either of the 
two vr^vr^vr'^ combinations have (a) a mass in the uj ot rj mass regions, and (b) all combinations for 
which the ■k^tt~tt^ invariant mass falls outside the uj and rj mass regions. The shaded region in (a) 
shows the corresponding distributions for the lo and rj sidebands (see text for signal and sideband 
region definitions). The superimposed curves are fits to a Gaussian signal function and a linear 
background as described in the text. The fitted numbers of events are shown at the top of each 
figure. 
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